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We investigate dynamic transitions and stochastic resonance phenomena in two-dimensional fully 
frustrated Josephson-junction arrays driven by staggered oscillating magnetic fields. As the temper- 
ature is lowered, the dynamic order parameter, defined to be the average staggered magnetization, is 
observed to acquire nonzero values. The resulting transition is found to belong to the same univer- 
sality as the equilibrium Z2 transition for small driving amplitudes while large driving fields appear 
to induce deviation from the universality class. The transition is also manifested by the stochastic 
resonance peak of the signal-to-noise ratio, which develops above the transition temperature. 

PACS numbers: 74.50. +r, 74.25.Nf, 05.40.-a 



The two-dimensional (2D) fully frustrated XY model, 
which provides a good description of the fully frustrated 
Josephson-junction array (FFJJA), has attracted much 
attention for the past decades, with regard to the Z2 sym- 
metry in .addition to the U(l) symmetry present in the 
system.Ert! Extensive studies have been devoted to the 
nature of the equilibrium transition in the system, lead- 
ing to general consensus on double transitions with some 
controversy as to the critical behavior of the Z2 transi- 
tioii.tm Recently, the kinetic Ising model, which is a pro- 
totype model with the Z2 symmetry, has been.-examined 
in the presence of oscillating magnetic fields-OB It has 
been revealed that spontaneous symmetry breaking takes 
place for finite strength of the oscillating fielder □ and that 
in two dimensions the transition belongs to the same unk 
versality class as the equilibrium 2D Ising transition.^ 
The signal-to-noise rajfcio (SNR) has also been found to 
exhibit double peaks J3 one above the transition temper- 
ature and the other below the transition temperature, 
the emergence of which has been explained in terms of 
the time-scale matching between the inherent relaxatio: 
time of the system and the external periodic driving 
Such interesting results in the kinetic Ising model invokes 
interest in the dynamics of the FFJJA, which possesses 
the same Z2 symmetry as the Ising model together with 
the additional continuous U(l) symmetry. Furthermore, 
unlike the Ising model, the FFJJA has real intrinsic dy- 
namics derived from the Josephson relations and thus 
grants direct experimental realizations^ this makes the 
study of the FFJJA more attractive. It is a challenge to 
prepare oscillating fields which couples relevantly to the 
Z2 degrees of freedom in the FFJJA. One possible real- 
ization may be provided by a periodic array of magnetic 
particlest2l which is placed .under the FFJJA and driven 
back and forth periodically.^ There have been some stud- 



ies on the dynamic properties of the FFJJA]I2Mi3 and 
very recently dynamic transitions were investigated in 
the FFJJA driven by uniform dc currcnt.lij However, no 
results are available yet about dynamic transitions and 
stochastic resonance (SR) phenomena in the presence of 
oscillating fields. 

This paper investigates the dynamic behavior of the 
FFJJA in the presence of the transverse magnetic field, 
which is staggered in space and periodic in time, with at- 
tention paid to the dynamic transition and SR phenom- 
ena. At low temperatures the chirality in the FFJJA dis- 
plays antiferromagnetic ordering, to which the staggered 
oscillating field plays the same role as the oscillating field 
in the kinetic Ising model. To describe the dynamic tran- 
sition to the antiferromagnetic ordering, we define the 
dynamic order parameter as the staggered magnetiza- 
tion averaged over one period of the (staggered) driving 
field, and examine its behavior with respect to the am- 
plitude and the frequency of the driving field as well as 
the temperature. For small driving amplitudes, the scal- 
ing plot strongly indicates that the transition belongs to 
the same universality as the equilibrium Z2 transition. 
On the other hand, as the driving amplitude is raised, 
apparent deviations can be observed in the scaling plot, 
suggesting a universality class different from the Z2 one. 
We also examine the SNR and find a peak above the 
transition temperature, disclosing the SR present in the 
system. The absence of the lower peak is discussed in 
relation to the relaxation time of the system. 

We begin with the equations of motion for the phases 
{<pi} of the superconducting order parameters in super- 
conducting grains forming an L x L square lattice. In 
the resistively-shunted-junctioiwpodel with the fluctuat- 
ing twist boundary conditions Ji-3 they read: 



1 



<Ydt 



sm 



■ A) + rjij 



= 0, 



(1) 



where the primed summation runs over the nearest neigh- 
bors of grain i. We have used the abbreviations <fnj = 



- A^j and r^j 



rj— Yj with Yi = (xi,yi) denot- 



ing the position of grain i, and expressed the energy 
and the time in units of hl c /2e and h/2eRI c , respec- 
tively, with the critical current I c and the shunt resis- 
tance R. The thermal noise current r/ij is assumed to be 
white satisfying (r7*i (*-Hr)»7fcJ (*)) = ^TS( T )^ikSji-SuSj k ) 
at temperature T and the dynamics of the twist variables 
A = (A x , A y ) is governed by the equations 



dAa 

dt 



-y 

L 2 ^ 



(2) 



where L is the system size, ^2uj\ denotes the summation 
over all nearest-neighboring pairs in the a-direction (a = 
x,y), and r) Aa satisfies (r)A a (t+T)r) Aa (t)) = (2T/L 2 )S(t). 
The gauge field Aij (= — Aji), which incorporates the ef- 
fects of the transverse magnetic field, takes the form 



for r v -\- x 

nxi + TT(~l) x '+y- f sin fit for r] = r- + y' (3) 

which corresponds to the combined uniform constant and 
staggered oscillating fields. The frustration, i.e., the flux 
in units of the flux quantum through the plaquette at 
the dual lattice site = r< + (1/2) (x + y) thus consists 
of the uniform full dc component and the staggered ac 
component of amplitude fo and frequency fi: 



fo = 5 + (-l) x<4 *7osinnt. 



(4) 



For the study of the transition associated with the Z2 
symmetry in the FFJJA, it is convenient to consider the 
chirality 



C(R, t) = sgn 



. p 



sin 



and the staggered magnetization 

m(t) = ^5^(-l)-+*'C(R,t), 



(5) 



(6) 



R 



where J^p denotes the directional plaquette summation 
of links around dual lattice site R. To probe the dy- 
namic transition in the presence of an oscillating field, 
we define the dynamic order parameter as the staggered 
magnetization averaged over the period of the fieldu 



Q 



fi 
2^ 



i(t)dt 



(7) 



In numerical simulations, the sets of equations of motion 
(PJ) and (0) have been integrated via the modified Euler 



method with time steps At — 0.05. We have followed an 
annealing schedule with the equilibration time at least 
500 periods at each temperature. 

In the fully frustrated array, it is well known that its Z2 
symmetry produces two kinds of antiferromagnetic chi- 
rality ordering at zero temperature. The staggered driv- 
ing field in Eq. (|4|), which favors one of the two ground 
states in the first half of a period and the other in the 
latter half, is expected to induce oscillations between the 
two ground states when the driving amplitude is suffi- 
ciently large. Figure [l] displays the time evolution of the 
staggered magnetization m(t) at zero temperature, with 
the initial condition m(t=0) = 1, for various amplitudes 
of the driving field. For /q < f c (which is equal to 0.5 re- 
gardless of the driving frequency), the system remains in 
the state with m = 1 even in the presence of the driving 
field. As the amplitude is increased beyond the critical 
value f c , the system is driven out of the m = 1 state, 
giving rise to oscillations in the staggered magnetization 
between m = ±1. It is interesting that the residence time 
in one state is different from that in the other. Accord- 
ingly, for fo > f c , although the dynamic order parameter 
is not unity, it still does not vanish. ^From Eq. (^), we 
expect the array to stay at the m = — 1 state only during 
the interval given by /osinfii < —1/2; this leads to the 
dynamic order parameter for fo > f c , 



0(/o) = l--cos- 1 -^, 

7T 2/o 



(8) 



which depends on the driving amplitude fo but not on 
the driving frequency fi. In Fig. g we plot the dynamic 
order parameter Q versus the amplitude fo at zero tem- 
perature. The numerical data indeed agree very well with 
the analytical results given by Eq. (P), establishing the 
validity of the analytical argument. 

At sufficiently high temperatures, the influence of the 
oscillating field may be neglected and the staggered mag- 
netization fluctuates randomly with time, resulting in the 
null value of the dynamic order parameter. Accordingly, 
we expect that there is a transition between the ordered 
state at low temperatures and the disordered state at 
high temperatures. Figure |3|, in which the behavior of 
(Q), the ensemble average of the dynamic order param- 
eter, with the temperature is plotted for various driving 
frequencies, demonstrates clearly the existence of such a 
transition: The dynamic order parameter, starting from 
zero at high temperatures, develops as the temperature 
is reduced. It then grows rapidly around a certain tem- 
perature, from which the transition temperature T c can 
be inferred, and saturates eventually as the temperature 
approaches zero. The saturation value assumed by the 
dynamic order parameter at zero temperature is unity for 
fo <lc as shown in Fig. |(a); for f > fc, it is given by 
Eq. (g) [see Fig. ||(b)]. It is indicated that in both cases 
the transition temperature becomes higher with the in- 
crease of the driving frequency fi. 
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The transition temperature in the thermodynamic 
limit can be convenientljj-obtained from the considera- 
tion of Binder's cumulantll3 



U L = 1- 



{Q 4 



3<g 2 ) 2 ' 



(9) 



This quantity is known to be independent of the system 
size at the transition temperature T c , which allows one 
to estimate the transition temperature by the crossing 
point of the cumulants for different sizes. [See the inset 
of Fig. in which the cumulant is plotted as a function 
of the temperature T for various system sizes.] In this 
manner we have estimated the transition temperature for 
various values of the driving amplitude fo and frequency 
£7, and show in Fig. |] the resulting phase diagram for 
the dynamic transition on the T—fo plane. For given 
frequency f2, the transition temperature decreases mono- 
tonically with the driving amplitude fo, similarly to the 
case of the kinetic Ising model. As is increased, on the 
other hand, the transition temperature also increases, ex- 
panding the ordered region on the T—fo plane. It is thus 
concluded that high-frequency driving helps to establish 
order whereas large-amplitude one tends to suppress or- 
der. It should be noted here that the cooling and heating 
curves for the dynamic order parameter do not exhibit 
any hysteresis in the strong-field regime. This appar- 
ently indicates the absence of a tricritical point, above 
which the system undergoes a discontinuous transition, 
and contrasts sharply with the kinetic Ising models un- 
der oscillating fields both in the mean-field versioru and 
in two dimensions B 

Regarding the nature of the transition, it is revealing 
to recall that the system at zero temperature, residing in 
one of the two (m = ±1) states for fo < f c , displays oscil- 
lations between the two states as fo is raised beyond f c . 
This suggests that the critical behavior near the transi- 
tion may change qualitatively at fo = f c . To probe such 
possibility, we consider the scaling law for the dynamic 
order parametertj 



(Q) = L-V v F((T-T c )L l l v ) 



(10) 



with appropriate scaling function F(x) and plot (Q)L I3 '' / 
versus (T — T C )L 1 / U in Fig. || for various sizes. Fig- 
ure ||(a) shows the best collapse of the numerical data for 
fo = 0.3, yielding the critical exponents (3 = 0.10 ± 0.02 
and v = 0.82 ± 0.05; these values provide good collapse 
in the whole range of fo < f c . Note the good agree- 
ment of the obtained exponents with those of the equi- 
librium transition^, which strongly suggests that the dy- 
namic transition at weak driving fields belongs to the 
same universality class as the equilibrium Z2 transition 
in the fully frustrated XY model. Similar conclusion 
was also reached in the kinetic Ising model under oscil- 
lating fields.lj For / > f c , on the other hand, the scaling 
plot with the same critical exponents, apparently does 
not collapse into a single curve. Such dispersion, which 
begins to show up around fo = f c , indicates deviation 



of the critical behavior from that of the equilibrium Z2 
transition, manifesting the expected change of the na- 
ture of the transition at fo = fc- Instead the data fit 
well to the scaling form in Eq. @ with j3 = 0.36 ± 0.08 
and v = 0.9 ± 0.1, as shown in Fig. |(b) for f = 0.8. 
These values of the exponents, estimated from the best 
collapse for fo > f c , turn out to be independent of fo 
within numerical errors. 

Finally, to explore the possibility of the SR phenomena 
in the system, we compute the SNR defined to be 



SNR= 10 log! 



(11) 



Here the signal S is given by the peak value of the power 
spectrum of m(t) at the driving frequency while the 
background noise level N is estimated by the average 
power spectrum around the signal peak. Figure ^ dis- 
plays the SNR versus the temperature (a) for weak driv- 
ing (fo = 0.1) and (b) for strong driving (/ = 0.6). For 
fo < fc shown in (a), the SNR peak emerges above the 
transition temperature and signals the presence of the SR 
in the system, which is reminiscent of the SR behavior 
of the kinetic Ising field in weak oscillating fields .0 For 
fo > fc one may expect monotonic decrease of the SNR 
with the temperature since the array has inherent oscil- 
lations of the staggered magnetization m(t) even at zero 
temperature. This expectation appears to be valid at 
low temperatures as shown in Fig. |6j(b) . As the temper- 
ature is raised, however, the SNR changes its behavior: 
Sharp decrease of the SNR turns into gradual increase 
seemingly at the transition temperature and into gradual 
decrease again above the transition temperature, giving 
rise to a small broad peak. 

Here it is of interest that the system exhibits no peak 
below the dynamic transition, which seems to be in con- 
tradiction to the prediction of the double SR peaks, one 
below and the other above the dynamic transition tem- 
perature.u The origin of the double peaks is the time-scale 
matching between the external periodic driving and the 
relaxation time associated with the transition in the sys- 
tem. To disclose the origin of the discrepancy, we ex- 
amine the relaxation of the system from random initial 
configurations at given temperatures, and estimate the 
relaxation time t by fitting Q(t) — (Q) to the exponen- 
tial form ~ exp[—i/T] in the long-time scale. The inverse 
of the resulting relaxation time for fo = 0.1 is shown in 
Fig. ^. As expected, the inverse relaxation time increases 
with the temperature above T c , which is consistent with 
the critical behavior 



(T - T c y 



(12) 



with the critical exponents given by the equilibrium val- 
ues v = 0.813 ± 0.005 and z = 2.17 ± 0.04.B Thus sup- 
ported again is the proposition that the dynamic transi- 
tion under weak driving fields belongs to the same uni- 
versality class as the equilibrium transition. At low tem- 
peratures below T c , on the other hand, it is observed 
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that the relaxation time is enhanced enormously. Due to 
such enhancement of the relaxation time below T c , the 
system fails to match the external time scale, which is 
estimated to be r _1 w 0.032 from the SNR peak above 
T c . The origin of the enhancement at low temperatures is 
not clear yet, but it can be attributed to the existence of 
a large number of metastable states associated with, e.g., 
domain- wall excitations, which cannot be observed in the 
mean-field system considered in Ref. |[ Indeed faceted 
domain walls have been pointed-put to cause slow growth 
in the FFJJA below Xjj w 0.3 JiH where the enhancement 
of the relaxation-time data is conspicuous in Fig. g This 
strongly supports our argument concerning the failure of 
time-scale matching. 

In summary, we have studied dynamic properties of 
the 2D FFJJA driven by staggered oscillating magnetic 
fields, paying attention to the dynamic transition and SR 
phenomena. To describe the dynamic transition to the 
antifcrromagnetic ordering of the chirality in the system, 
we have introduced the dynamic order parameter, which 
is the staggered magnetization averaged over one period 
of the staggered driving field, and examined its behavior 
with respect to the amplitude and the frequency of the 
driving field as well as the temperature. The transition 
temperature has been computed by means of Binder's cu- 
mulants for various driving amplitudes and frequencies, 
and as a result, the phase diagram displaying boundaries 
of the dynamic transition has been obtained on the plane 
of the temperature and the driving amplitude for various 
driving frequencies. While the scaling analysis for weak 
driving has demonstrated the universality of the equi- 
librium Z2 transition, a universality class different from 
the Z 2 one has been suggested for strong driving. We 
have also examined the SNR and found a peak above 
the transition temperature, disclosing the SR present in 
the system. The absence of the lower peak has been dis- 
cussed in relation to the relaxation time of the system 
and attributed to the striking increase at low tempera- 
tures. 
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FIG. 1. Time evolution of the staggered magnetization at 
zero temperature for driving frequency Q./2n — 0.02 and am- 
plitudes fo — 0.4,0.6, and 1.2 from above. For clarity, the 
data corresponding to fo = 0.4 and 1.2 have been shifted by 
1.5 upward and downward, respectively. 
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FIG. 2. Dynamic order parameter as a function of /o at 
zero temperature in the system of size L = 16 and with fre- 
quency fl/2-K = 0.02. The dashed line represents the analytic 
result in Eq. (|). 
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FIG. 4. Dynamic phase diagram on the T—fo plane for var- 
ious driving frequencies D./2n = 0.02(D), 0.04(O), 0.08(A). 
The boundaries are determined by the crossing points of 
Binder's cumulant for several sizes L — 8, 16, 24, and 32. In- 
set: Binder's cumulant Ul as a function of the temperature 
for fo = 0.3 and Q./2n — 0.08 in the system of size L = 8(D), 
16(0), 24(A), and 32( V ). 
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FIG. 3. Dynamic order parameter as a function of temper- 
ature T in the system of size L = 16 for various driving fre- 
quencies n/2-K = 0.02(D), 0.04(0), 0.08(A) and (a) f = 0.3; 
(b) f = 0.6 . 
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FIG. 5. Scaling plot of the dynamic order parameter ver- 
sus the temperature for size L — 8(D), 16(0)> 24(A), and 
32(v) with Q/2n = 0.08 and (a) fo = 0.3; (b) fo = 0.8. The 
best collapse is obtained with (a) f3 = 0.10 and v — 0.82; (b) 
13 = 0.36 and v = 0.90. 



(b) 

FIG. 6. Signal-to-noise ratio versus the temperature in the 
system of size L = 16 for (a) fo = 0.1 and 0/27T = 0.02; (b) 
fo — 0.6 and Q/2tt — 0.08. The arrows indicate the positions 
of the transition temperatures. 
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